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Abstract 

By employing the QCD factorization approach, we calculate the new physics contributions 
to the branching radios of the two-body charmless B — ► PP decays in the framework of the 
minimal supergravity (mSUGRA) model. Within the considered parameter space, we find that 
(a) the supersymmetric (SUSY) corrections to the Wilson coefficients C/% (k = 3 — 6) are very 
small and can be neglected safely, but the leading order SUSY contributions to CV 7 (Mw) and 
Cs g (M\v) can be rather large and even change the sign of the corresponding coefficients in 
the standard model; (b) the possible SUSY contributions to those penguin-dominated decays 
in mSUGRA model can be as large as 30 — 50%; (c) for the well measured B — ► Kir decays, 
the significant SUSY contributions play an important rule to improve the consistency of the 
theoretical predictions with the data; (d) for B — > Kn' decays, the theoretical predictions of the 
corresponding branching ratios become consistent with the data within one standard deviation 
after the inclusion of the large SUSY contributions in the mSUGRA model. 
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I. INTRODUCTION 



As is well known, the precision measurements of the B meson system can provide an 
insight into very high energy scales via the indirect loop effects of the new physics beyond 
the standard model (SM) [H0]. Although currently available data agree well with the SM 
predictions, we generally believe that the B-factories can at least detect the first signals 
of new physics if it is there. 

Among the B — ► PP (P stands for the pseudo-scalar light mesons) decay channels 
considered in this paper, twelve of them have been measured with good accuracy. And 
the data indeed show some deviations from the SM expectations: 

• The Kr)' puzzle, the observed B — > Kr( branching ratios 0, 0, HI are much larger 
than the corresponding SM predictions, appeared several years ago, and there is 
still no convincing theoretical interpretation for this puzzle after intensive studies 
in the framework of SM [f| and the new physics models^]. 

• The Ktt puzzle comes from the ratios R c and R n for the four well measured B — > Kit 
decay rates as defined in Ref. The SM prediction is R c = R n by neglecting the 
small exchange- and annihilation-type amplitudes [H, while the present data P 
yields 



2T(B + -> K~ 



7T 



0\ 



R c XP = T#^^f = 1-15 ±0.13, (1) 

2r ( ^ W r a78±al0 (2 > 

A discrepancy of 2.8a exist here. 

• For B — > 7T°7r decay, the measured branching ratio Br(B — ► = (1.9 ± 0.5) x 

io" 6 11111 is about five times larger than the SM prediction. 

Although not convincing, these discrepancies together with the so-called <j >K K anomaly 
(Toj may be the first hints of new physics beyond the SM in B experiments (ill . 

Up to now, the possible new physics contributions to rare B meson decays have been 
studied extensively, for example, in the Technicolor mode ls|l3lj . the two-Higgs-doublet 
modelsjli], [l5| and the supersymmetric models [3, 17, EM El- Among the various new 



physics models, the supersymmetric models are indeed the most frequently studied models 
in searching for new physics in B meson system. The minimal supersymmetric standard 
model (MSSM)[2(j is the general and most economical low-energy supersymmetric exten- 
sion of the SM. But it is hardly to make definite predictions for the physical observables 
in B meson decays since there are more than one hundred free parameters appeared in the 
MSSM. In order to find the possible signals or hints of new physics beyond the SM from 
the data, various scenarios of the MSSM are proposed by imposing different constraints 
on it j2£|. The minimal supergravity (mSUGRA) model (2l| seems to be a very simple 
constrained MSSM model, since it have only five free parameters tan/3, mi , mo, Ao, and 
sign(n) at the high energy scale. 

The previous works in the framework of mSUGRA model focused on the semileptonic, 
leptonic and radiative rare B decays. In Refs. fl6l 17, [3, 22||, for example, the authors 
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studied the rare decays B — > X s ^, — > X s ll, B — > Z + / and the 5° — 5° mixing in the 
mSUGRA model, and found some constraints on the parameter space of this model. 

For B — > PP decays, they have been studied in the SM El 0, H H B 13 , the 
Technicolor models[13] and the two-Higgs-doublet models[14|]. In Ref.|3]J, Mishima and 
Sanda calculated the supersymmetric effects on B — > <fik decays in the PQCD approach j30l] 
and predicted the values of CP asymmetries with the inclusion of the supersymmetric con- 
tribution. In this paper, we calculate the supersymmetric contributions to the branching 
radios of the twenty one B — > PP decay modes in the mSUGRA model by employing 
the QCD factorization approach (QCD FA) jjjl. I2H l2?|. The contributions from chirally 
enhanced power corrections and weak annihilations are also taken into account. We find 
that the branching ratios of some decay modes can be enhanced significantly, and these 
new contributions can help us to give a new physics interpretation for the so-called u Krj' n 
puzzle. 

This paper is organized as follows. In section II, we give a brief review for the minimal 
supergravity model. In section III, we calculate the new penguin diagrams induced by 
new particles and extract out the new physics parts of the Wilson coefficients in the 
mSUGRA model. The calculation of B — > PP decays in QCD factorization approach 
is also discussed in this section. In section IV, we present the numerical results of the 
branching ratios for the twenty one B — > PP decay modes in the SM and the mSUGRA 
model, and make phenomenological analysis for those well measured decay modes. The 
final section is the summary. 



II. OUTLINE OF THE MSUGRA MODEL 



In the MSSM, the most general superpotential compatible with gauge invariance, renor- 
malizability and R-parity conserving is written as 2ol |: 



'-a/3 



(3) 



where fz>, fu and f E are Yukawa coupling constants for down-type, up- type quarks, 
and leptons, respectively. The suffixes a, (3 = 1,2 are SU(2) indices and i,j=l,2,3 are 
generation indices, e a p is the antisymmetric tensor with e 12 = 1. In addition to the SUSY 
invariant terms, a set of terms which explicitly but softly break SUSY should be added 
to the supersymmetric Lagrangian. A general form of the soft SUSY-breaking terms is 
given as 20]: 

- C sof t = (m 2 Q ).. qtttLj + {ml).. u Ri u Rj + (m 2 D ) .. d* Ri d Rj + (m 2 L ).. l+f Lj 
+ (m 2 E ).. e* Ri e Rj + A 2 h+hi + A 2 2 h+h 2 



A uijq2Mu Rj + A Dij hff Li d* Rj + A Eij h*P Li e Rj + B^h p 2 



+ ^m B BB + ^m w WW + ^m d GG + H.C. 



(4) 



where qu-, u* Ri , d* Ri , In, e* Ri , and hi and h 2 are scalar components of chiral superfields 
Qi, Ui, Di, Li, Ei, Hi, and H 2 respectively, and B, W, and G are U{1)y, SU{2)l, and 
SU(3)c gauge fermions. And the terms appeared in Eq.(@J) are the mass terms for the 
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scalar fermions, mass and bilinear terms for the Higgs bosons, trilinear coupling terms 
between sfermions and Higgs bosons, and mass terms for the gluinos, Winos and binos, 
respectively. 

In the mSUGRA model, a set of assumptions are added to the MSSM. One underlying 
assumpsion is that SUSY-breaking occurs in a hidden sector which communicates with the 
visible sector only through gravitational interactions. The free parameters in the MSSM 
are assumed to obey a set of boundary conditions at the Plank or GUT scale: 



«1 = «2 = «3 = OCX, 

(™q)« = ( m u)ij = {m 2 D )ij = (mDij = (m|) y = (mg)% 
A 2 = A 2 = m 2 , 

AlJij — fllijAo, Aoij = foijAo, Afiij = f'EijAo, 

nig = = m^ = mi (5) 

where a« = g 2 /(47r), and gi (i=l,2,3) denotes the coupling constant of the U(1)y, SU{2)l, 
SU(3)c gauge group, respectivel y. The unification of them is verified according to the 
experimental results from LEPl|32j and can be fixed at the Grand Unification Scale 
Mq UT ~ 2 x 10 16 Gef . Besides the three parameters mi, m and A , the super symmetric 
sector is described at GUT scale by the bilinear coupling B and the supersymmetric 
Higgs(ino) mass parameter /i. However, one has to require the radiative electroweak 
symmetry-breaking (EWSB) takes place at the low energy scale. The effective potential 
of neutral Higgs fields at the tree-level is given by (to be precise, one-loop corrections to 
the scalar potential have been included in the program we used later) 



Vm gg s = m\\h\\ 2 + m\\hl\ 2 + ml{h\hl + H.C) 



where we have used the usual short-hand notation: m\ = (/i 2 + Af), m\ = (fi 2 + A 2 ,), 
m\ = B^i. The radiative EWSB condition is 

,9V, , dV . 

where the value h® denotes the vacuum expectation values of the two neutral Higgs 
fields as = v cos/3, (h®) = wsin/3 with v = \7AGev. From Eq.((7j), we can determine 
the values of fi 2 and Bfi: 

fi 2 = i [tan 2/? (A 2 , tan /3- A 2 cot /3) - M 2 Z \ 

Bfx = ^sin2/?[A 2 + A 2 + 2^ 2 ] (8) 

Through Eq.fjHJ) we can see the sign of \x is not determined. Therefore only four continuous 
free parameters, and an unknown sign is left in the mSUGRA model. They are: 

tan mi , mo, A Q , sign(fi) (9) 
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In the mSUGRA model, all other parameters at the electroweak scale are then deter- 
mined through the five free parameters by the GUT universality and the renormalization 
group equation (RGE) evolution. In this paper, we calculate the SUSY and Higgs particle 
spectrum through a Fortran code: SUSPECT version 2.1 |33|. The important features 
of this code include (a) the renormalization group evolution between low and high en- 
ergy scales; (b) consistent implementation of radiative electroweak symmetry breaking; 
and (c) calculation of the physical particle masses with radiative corrections. Using this 
code, we obtain the SUSY and Higgs particle masses, and the mixing angles of squarks 
at the electroweak scale. From these Low-energy supersymmetric parameters, the mixing 
matrices T u ,T D for the up-type and the down-type squarks, the mixing matrices U, V, N 
for charginos and neutralinos are determined. The explicit expressions of the two 6x6 
mixing matrices T u and T D , two 2x2 matrices U and V, and a 4 x 4 matrix N can be 
found in Refs. 0, . 



III. THE BASIC THEORETICAL FRAMEWORK FOR B PP 

In this section, we present the theoretical framework and the relevant formulas for 
calculating the exclusive nonleptonic decays of the B ± and B° mesons into two light 
pseudoscalar mesons. 



A. Effective Hamiltonian and relevant Wilson coefficients in SM 

In the SM, if we take into account only the operators up to dimensions 6, and assume 
m& ^> m s , the effective Hamiltonian for the quark level three-body decay b — > qq q 
(q G {d, s}, q G {u, d, s}) at the scale \i reads [36( 

C { 2 

10 >| 

~V tb V* ]T 0^)0^) - V tb V* [C 7 M0 7 M + C 8g ^)0 8g ^)} (10) 

j=3 J 

where V^V^ is the products of elements of the Cabbibo-Kabayashi-Maskawa quark mixing 
matrix [3 71] . And the current- current (Oi )2 ), QCD penguin (0 3: 4^ :6 ), electroweak penguin 
(0 7>8 ,9,io), electromagnetic and chromomagnetic dipole operators (0 77 and 8g ) can be 
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written as 



Ox = (qu)v-A(ub) V -A, 0% = (q a Up)v-A(upb a ) v ^ A , 

0{ = (qc) V - A {cb)v-A, 0% = (q a cp) V - A {cpb a )v-A, 

3 = {qb) V -A ^2(qq')v-A, A = (q a bp) v „ A ^(qp q' a )v-A, 

<7 Q 



5 = (qb) v ^ A ^2(q'q')v+A, 6 = (q a b(s)v-A ^(qp q' c 



a)V+A, 



1 1 

3, ^ , , , x „ 3 



O7 = 7;(qb)v-A^e q >(q'q') v+A , 8 = -{q a bp) V -A^e q '{qp 'q' a )v+A, 



q 



O9 = -{qb) V - A ^e q > (q q')v-A, O 10 = ^{q a bp) V - A ^ e q > (qp'q' a )v-A, 

q q 

77 = ^m b q a a^(l + l5 )b a F^ 8g = JLm^" (1 + 75 )T^G%, (11) 

where T a (a = 1, 8) stands for SU (3) c generators, a and (3 are the SU (3)c color indices, 
and V ± A = 7^(1 ± 75) by definition. The sum over q' runs over the quark fields that 
are active at the scale \i = 0(mt,), i.e., q' G {u, d, s, c, b}. 

To calculate the non-leptonic B meson decays at next-to-leading order in a s and to 
leading power in AQCD/mb, we should determinate the Wilson coefficient Ci(Mw) through 
matching of the full theory onto the five-quark low energy effective theory where the W ± 
gauge boson, top quark and all SUSY particles heavier than M w are integrated out, and 
run the Wilson coefficients down to the low energy scale /i ~ 0{m^) by using the QCD 
renormalization group equations. In table |1J we simply present the numerical results of 
the LO and NLO Wilson coefficient in the NDR scheme in different scales. More detailed 
analytical expressions can be found for example in Refs. 3^, 38| . 



B. Wilson coefficients in the mSUGRA model 

In the mSUGRA model, the new physics contributions to the rare decays will manifest 
themselves through two channels. One is the new contributions to the Wilson coefficients 
of the same operators involved in the SM calculation, the other is to the Wilson coefficients 
of the new operators such as operators with opposite chiralities. In the SM, the latter is 
absent because they are suppressed by the ratio m s /mb- In the mSUGRA model, they 
can also be neglected, as shown in Ref . 39] . Therefore we here use the same operator base 
as in the SM. 

It is well known that there is no SUSY contributions to the Wilson coefficients at the 
tree level. There are five kinds of contributions to the quark level decay process b — > qqq 
at one-loop level, depending on specific particles propagated in the loops: 

• the gauge boson W ± and up-type quarks u, c, t, which leads to the contributions in 
the SM; 

• the charged Higgs boson H ± and up-type quarks u, c, £; 
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TABLE I: In the NDR scheme, the values of LO and NLO Wilson coefficients Ci(n) for 
\i = mb/2,mb,2mb- Input parameters being used are A^I = 0. 225Gev, sin 2 9\y = 0.23, nib = 
4.62Gev,m t = 175Gev,M w = 80AGev, and a em = 1/128. 





H = m b /2 


/* = 


m b 


V = 


2mb 




LO 


NLO 


LO 


NLO 


LO 


NLO 


c x 


1.179 


1.134 


1.115 


1.080 


1.072 


1.043 


c 2 


-0.370 


-0.280 


-0.255 


-0.180 


-0.171 


-0.104 


c 3 


0.019 


0.020 


0.012 


0.013 


0.008 


0.008 


Ca 


-0.037 


-0.048 


-0.027 


-0.034 


-0.018 


-0.023 


C 5 


0.010 


0.012 


0.008 


0.010 


0.006 


0.007 


c 6 


-0.050 


-0.062 


-0.033 


-0.040 


-0.021 


-0.026 




0.018 


-0.008 


0.028 


0.007 


0.046 


0.030 


C%ja em 


0.055 


0.055 


0.035 


0.035 


0.023 


0.023 


C%ja em 


-1.398 


-1.420 


-1.318 


-1.337 


-1.255 


-1.270 


Cio/a em 


0.415 


0.395 


0.286 


0.273 


0.191 


0.183 


C77 


-0.360 


-0.334 


-0.316 


-0.307 


-0.281 


-0.282 


C 8 g 


-0.167 




-0.150 




-0.136 





• the charginos xf 2 an d the scalar up-type quarks u, c, t; 

• the neutralinos X? 2 3 4 an d the down-type squarks d, s, b; 

• the gauginos g and the down-type squarks d, s, b. 

The new physics contributions from those superparticle loops may induce too large 
flavor changing neutral currents (FCNCs). To escape from the so-called SUSY flavor 
problem, degeneracy of masses of squarks and sleptons among different generations has 
been assumed in the minimal SUGRA model. 

In order to determine the new physics contributions to Wilson coefficients Ci(i = 
3,4,5,6), CV 7 , and C& g (we ignore the new physics contributions to Ciii = 7,8,9,10) 
because they are suppressed by a factor of a em /a s ) at the M w scale, we need to calculate 
the Feynman diagrams appeared in Fig^ First, by employing conservation of the gluonic 
current, we can define the effective vertex of the b — >• qg penguin processes as in Ref . |4oj : 

r> 2 ) = ^u q (p q )T a V,(q 2 )Mp b ) (12) 

with 

VM 2 ) = (q 2 9,u-q,qu)Y[FiL(q 2 )PL + F 1R (q 2 )P R \ 

+i<y^q" [F 2L (q 2 )PL + F 2R (q 2 )P R ] (13) 

where F\(q 2 ) and F 2 {q 2 ) are the electric and magnetic form factors, q = pb — p q is the 
gluon momentum, and Pl(r) = (1 ±75)/2 are the chirality projection operators. 
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By calculating the Feynman diagrams as shown in FigsJ^b)-(f), we find (in the naive 
dimensional regularization (NDR) scheme,) the new physics parts of the Wilson coeffi- 
cients at the scale M w 

-i -l 



cr(Mw) 

Cg P {M w ) 



a a (Af, 



w 



24tt 
F 2 r(0) 



G] 



■X 



G p 



y/2 



V2 
X t m b 



A k F 1L (0) 



-1 



(14) 
(15) 



where A k = {— 1, 3, — f , 3} for k = {3, 4, 5, 6}, and \ t = V t * q V t b- In addition, since q 2 <C m~ 
where m q is the mass of the heavy scalar fermions, we set q 2 = for the form factors F\$ 
as an approximation l . The explicit expressions of the form factors F 1L (0) and F 2 r(0) 
induced by supersymmetric particles are the following 



*S~(o) 

*£"(0) 

f!m 

*?*(0) 



G p 

'F 

V2 



\ t xth cot 2 (3f 5 (x t h), 



Gp 

—=\ t m b x th [cot 2 Pfi(x th ) + f 3 {x th )] , 



9l 



4m 



(16) 
(17) 

(18) 



9 1=1 

2 6 
9 s 



d \ 3 



1=1 



d \3 



G 



2 6 



EE^U! ^Cl) 1 ^ ( T CL)T ft ( X x-uj) 
VZ a=l 1=1 

2 6 

V2G F J2J2 X ^ ( r CLY aj \FcL)Tm*M 



(19) 
(20) 



a=l 1=1 
^d \ a3 



X a Uj) 



+ ( t cr) i m x-h( X x-uj) 



G 4 6 

^lL( ) = 1= ^ y^^H"dr( r ^)a?( r iVL)j 3 /7(^0d r ) ; 
v q=1 /=l 



(21) 
(22) 



4 6 



(0) = -V^Gjp^E^W^ ( r ^)f^/2(^ 



a=l 1=1 

■^d \ce3„ 



+(r% R )rm,oj 4 (x- OJi : 



(23) 



where j = 1 for 6 — > d and j = 2 for 6 — > s decay, respectively. Xjj = mf/m 2 and m« 
is the mass of the particle i. In our calculations, we have set m q = for g = u,d,s 



1 Eos. (|f 4(1 and (|15|l differ from those appeared in Ref.|41j, but our final analytic expressions for 



C£L P (Mw) are the same as that in Ref.|l8l| except for the definition for Wilson coefficients. 
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w~ 

b 1 u,c,t x s, d 
». I » , > I > — 



5 ' u,c,t s s, d 
» I — » a » 1 » 



ff, 7 (6) 



6 ' 9 9 \ s, d 
» ' — » ) > I — 



(c) 





.9, 7 ( e ) 




ff- 7 (/) 



FIG. 1: Five kinds of loop contributions to 6 — ► qg,qj- (a) SM contribution, (b) charged Higgs 
contribution ; (c-d) gluino contribution; (e) chargino contribution; (f) neutralino contribution. 



since mj 3> m Ul md and m s . The one-loop integration functions fi(x) and the coupling 
constants Fq( L r \, ^c(l r)> ^n(l r) which appear in F iL and F 2 r are listed in Appendix 1X1 
Using the form factors in Eqs. (jl(ij) - (j23j) . we obtain the analytic expressions for C^ p (Mw) 
and Ca 3 P (Mw). The Wilson coefficient C^ g p (Mw) as given in EchflTK]) is the same as 
that in Ref. [18[ except for some differences in expression. In Ref.[l8|, the CKM factor 
— Xt has not been extracted from Wilson coefficients, and the CKM matrix elements have 
been absorbed into the definition of the coupling constant 1^. See Appendix El for more 
details. 

For the effective vertex of the supersymmetric b — > (77 penguin processes, we only 
consider it's contributions to CV 7 . The explicit analytical expressions of the SUSY con- 
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tribution to CV 7 induced by new particles have been given in Ref. [1 



C 7 H ;(m w ) 
(m w ) 



~^ x th { cot 2 /3 



;fi{x t h) + f2(x th ) 



+ 



\h{xth) + U{x t h) 



(r^);/ 2 (^) + (r^) 7 ^/4fe; 

III'}) 



(24) 
(25) 



a=l 1=1 



a3 



(26) 



c? 7 W) = -^EE%*( r K| r u; 3 / 2 (^,) 



Q=l 7=1 



+ ( r 7V J R)/ 3 — 7^/4(^0 



(27) 



Now, we found all the supersymmetric contributions to the relevant Wilson coefficients. 
We should remember that, the only source of flavor violation in the mSUGRA model is 
the usual CKM matrix in the SM. The flavor violation in the sfermion sector at the 
electroweak scale is generated radiatively in the mSUGRA model and consequently small. 
Therefore, If we take the mixing matrices T u and T D as given in the Appendix of Ref.[34T| 





1 



















1 






















f 



















t 






















cos 9f 








sin B~ t 










cos 6^ 








sin Q~ b 













f 



















f 






















f 



















f 













— sin Oi 








cos 61 










— sin 9i 








cos d~ h 



(2f 



the gluino- and neutralino-mediated diagrams will not contribute to the decay processes 
considered here. The new physics contributions will come from the charged-Higgs and 
chargino diagrams only. 



C. B — > PP decays in QCD factorization 

To calculate the decay amplitude of the processes B — > PP, the last but most impor- 
tant step is to calculate hadronic matrix elements for the hadronization of the final-state 
quarks into particular final states. At the present time, many approaches have been put 
forward to settle the intractable problem. Such as the native factorization 42l |. the gen- 



eralized factorization HHJ, the QCD FA 0, and the PQCD approach [ij. In 
this paper, we employ the QCD FA to calculate the branching ratios of B — ► PP decays. 

In QCD FA, the contribution of the non-perturbative sector is dominated in the form 
factors of B — > P transition and the nonfactorizable impact in the hadronic matrix 
elements is controlled by hard gluon exchange. In the heavy quark limit m& ^> Aqcd and 
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to leading power in AQCD/f^b, the hadronic matrix elements of the exclusive nonleptonic 
decays of the B meson into two light pseudoscalar mesons Pi, P 2 (Pi absorbs the spectator 
quark coming from the B meson) can be written as [25| 

(P 1 P 2 \O i \B) = V^ Pl / dxTl j (x)^ P2 (x) + (P 1 ^P 2 ) 

+ [ d£ / dx I dyT l II ^,x,y)<l> B (0 ( l>p 1 (^P 2 (y) (29) 
Jo Jo Jo 

where F^ Pl is the form factor describing B — > Pi decays. and T/ 7 denote the per- 
turbative short-distance interactions and can be calculated by the perturbation approach. 
<f>x( x ){X = B,Pi 2 ) are the universal and nonperturbative light-cone distribution ampli- 
tudes (LCDA) for B and Px t2 meson respectively 26]. Weak annihilation effects are not 
included in Eq. flUp . 

Consider the low energy effective Hamiltonian Eq.()10|) and the unitary relation of the 
CKM matrix, the decay amplitude can be written as 

A(B -> P 1 P 2 ) = % X) ^V^Ciifi^P^lB) (30) 

p=u,c i 

the effective hadronic matrix elements (PiP 2 \Oi\B) can be calculated by employing the 
QCD factorization formula Eq. (J29|) . When considering order a s corrections to the hard 
scattering kernels Ty and T/ 7 from nonfactorizable single gluon exchange vertex correction 
diagrams, penguin diagrams and hard spectator scattering diagrams and the contributions 
from the chirally enhanced power corrections 2 , Eq. (j30|) can be rewritten as|29j 

A*{B - P^) = % E EWXM^IO^}, (31) 

* p=u,c i 

Here {P\P 2 \Oi\B) p is the factorized matrix element and can be factorizes into a form factor 
times a decay constant. The explicit expressions for the decay amplitudes of B — > P\P 2 
can be found in Ref . [23l| . For the processes involving rf^ in the final states, Ali et al. 
[23l | included the terms directly proportional to the so-called charm decay constant f c ^ 

of the 7/') meson in the decay amplitudes. We here ignored these terms because they 
are very small in size. For the charmless B meson decays considered here, the hadronic 
matrix elements {PiP 2 \01 2 \B)p have no contributions. Following Beneke et al. |26|, every 
coefficient ai{Pi, P 2 )(i = 1 to 10) can be split into two parts: 

Oi(Pi, P 2 ) = a f)J (Pi, P 2 ) + cvKPx, P 2 ) (32) 



2 For more details of various contributions and the corresponding Feynman loops, see for example 
Refs. tifl l22l and references therein. 
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with 



n C 2 C F a s C 2 C F Tia s 

ay = Ci + — [1 + -^Vy , «y/ = ^^^V- 

a V = C* + ^[1 + - ^ = A^HVT^' 

CV CVa s C±C F ira s 

as,/ = C 3 + ^-[1 + - a 3 ,,/ = —-^-H^, 

a ^ ~ Ci + W c [1 + — Vp2] + a4 - 7/ " iv^nvr^ PlP2 ' 

a 5 , = c 5 + ^[i + -^ri-vpji *w = w c ^r { - H 

<i = ^ + jA 1 - 6-^-1 + -^x- = 0, 



C8 ri , C F a s . t / m C$ C F 7ra & 



a ' J = C " + IV? 1 + ~^ J ' ° w ' = ATHvT' fffi ' p " 

<o,, = C + |[! + °gv n] + ^5f, ^K^Hn, (33) 



where N c = 3, Cp = 4/3, a^j = di,i{p) and a^// = aiji{jj, h ) with /i ~ m fe and /i^ = \fKhjl 

with Aft = 0.5Get> as in Ref . . The terms V P ^ result from the vertex corrections, H ^ P P 

describe the hard-scattering spectator contributions, Pp 2 2 and P P2 3 {P P E ^ and P P E ^ ) 
arise from the QCD (electroweak) penguin contributions and the contributions from dipole 
operator 0$ g ( Oj-y ) . For the four penguin terms, the subscript 2 or 3 indicates the twist of 

the corresponding projection. The explicit expressions of the functions V P \ H P ^ P2 , P P2 2 , 

P P2 3 , P p E 2 W and P P2 E ^ can be found in Ref. jifj] . 

In QCD FA, the non-factorizable power-suppressed contributions are neglected. How- 
ever, the hard-scattering spectator interactions and annihilation diagrams cannot be ne- 
glected because of the chiral enhancement. Since they give rise to infrared endpoint singu- 
larities when computed perturbatively, the y ca n only be estimated in a model-dependent 
way and with a large uncertainty. In Refs. |26l |2tJ these contributions are parameterized 
by two complex quantities, Xh and Xa, 

X H . A =(l + p H . A e l ^)\n^ (34) 

Aft 

where Aft = 0.5 GeV, 4>h,a are free phases in the range [—180°, 180°] and Ph,a are real 
parameters varying within [0, 1]. In this paper, we use the formulas as given in Ref.jiil] 
directly to estimate the annihilation contributions to specific final state. Under the con- 
vention of Ref. |26[ , the annihilation amplitude can be written as 

A ann {B _ K <h £ £ V^/a/ft/ft^l, A) (35) 

p=u,c i 
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where fu are the decay constants of B meson and final-state hadrons respectively. 
The coefficients b{{P\,P-i) describe the annihilation contributions and generally depend 
on quantity X^. For explicit expressions of coefficients 6, one can see Ref . (26| . 
Now the total decay amplitudes can be written as 

A{B -> P 1 P 2 ) = A f {B -> P X P 2 ) + A ann {B -> PiP 2 ) , (36) 

the corresponding branching ratio then takes the form 

B{B -> PP 2 ) = r B ^tI^(5 - PiP 2 )| 2 , (37) 

where tb is the B meson lifetimes, and \P C \ is the absolute values of two final-state hadrons' 
momentum in the B rest frame. For the CP-conjugated decay modes, the branching ratios 
can be obtained by replacement of A p — > A* in the expressions of decay amplitudes. 

The new physics contributions to the branching ratios of B —>■ PP decays will be 
included by using the Wilson coefficients Cj with the inclusion of the new physics parts 
as described in Eos. (11411151) and (I24H27I) . 



IV. NUMERICAL CALCULATIONS 

In this section, we first give the input parameters needed in numerical calculations, 
and then present the numerical results and make some theoretical analysis. 



A. input parameters 



The parameters (A, A, p, fj) in Wolfenstein parametrization of the CKM matrix. At 
present, the parameter A and A have been well determined by experiments. In 
numerical calculation, we will use A = 0.854, A = 0.2196, p = 0.22 ± 0.10, and 
f] = 0.35 ± 0.05 as given in Ref.Q. 

Quark masses. When calculating the decay amplitudes, the pole and current quark 
masses will be used. For the former, we will use 

m u = 4.2Mev, m c = 1.5Gev, m t = 175Gev, 

m d = 7.6Mev } m s = 0.122Gew, m b = 4.Q2Gev. 



The current quark mass depends on the renormalization scale. In the MS scheme 
and at a scale of 2GeV, we fix 

m u {2Gev) = 2AMev, m d {2Gev) = QMev, 

m s {2Gev) = 105Mev, m b {mE) = 4.26Gev, 
as given in PDG 2002 and then employ the formulae in Ref. 



m(fj) = m(p ) 



a s (po) 



1M 



1 + 



Jm 

Wo 



Pi 7m j a s (p) - a s (p ) 



(38) 
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to obtain the current quark masses at any scale. The definitions of a s , 7^ , 7m\ Po, 
and Pi can be found in Ref.|38|. 



Form factors and decay constants. Following Ref.[27], we also use 

F<f ^(0) = 0.28 ± 0.05, F M (0) = 0.34 ± 0.05. (39) 

The decay constants of ir, k and B are |26| 

f n = UlMev, f k = 160Mev, f B = 180Mev 

For 77 and rj ', mixing happens between them. The decay constants of them can be 
parameterized by f q , f s and the mixing angle of 77 — r\ 

# = ft = f%cos<i>, f° = -f s sincf> 

with 

/, = (1.07±0.02)/ ff , / s = (1.34±0.06)/ w , = (39.3 ± 1.0)° 
Similarly, the form factors ^^^(0) and Fq~* v (0) are parameterized as in Ref . |43| . 

For the parameters pn,A and 4>h,Ai we do not consider the variation of these param- 
eters but fix 

Pa — 0.05, A = 1O°, p H = 0, <p H = 0° 

in numerical calculation. For the parameter A b appeared in the B meson light-cone 
distribution amplitude, we also take Xb = (350 ± 150) Mev as in Ref. |26|| . 



For the well-known tt, K, rf^ and B meson masses, as well as the B meson lifetimes, 
we use the values as given in Ref. [32| • 



• The SUSY parameters at electroweak scale. Within the parameter space still allowed 
by known constraints from the data (IR fl8L E3 ] ( such as the strong constraints from 
the precise measurements of BriB — > X 8 j)), we choose two sets of SUSY parameters 
of the mSGURA model at the high unification energy scale as listed in Table |H] 
The resulting masses of charged Higgs boson and charginos obtained by using the 
program SUSPECT V 2.1 [33j are also given in Table ITT1 

In numerical calculations, we always use the cental values of above input parameters 
unless explicitly stated otherwise. 



B. Wilson coefficients: Case A and B 

From explicit calculations, we find that the SUSY corrections to B — > PP decays are 
mostly induced by the new physics parts of the CV 7 and Cg p , while the coefficients C^ p 
(k = 3,4,5,6) are indeed too small to modify their SM counterparts effectively. The 
numerical results show that the C 7l (m b ) and C 8ff (mb) in mSUGRA model can be quite 
different from that in the SM, and can even have the opposite sign compared with their 
SM counterparts. 
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TABLE II: Two sets of SUSY parameters to be used in numerical calculation. And the corre- 
sponding mass spectrum of charged Higgs boson and the charginos. All masses are in unit of 
GeV. 



Cases 


tan (3 


mi 

2 


m 


A 


sign(fi) 


m H ± 




m xt 


Case- A 


2 


300 


300 







782.3 


247.0 


595.9 


Case-B 


40 


150 


369 


-400 


+ 


330.2 


109.6 


312.3 



1. Case A 



We firstly consider the Case A. For the SUSY part, since we take the mixing matrix T u 
and T D as given in Eq. (J28|) . the Feynman diagrams induced by the gluino and neutralino 
exchanges do not contribute to the quark level decays b — > (s, d)^ and b — > (s, d)g. To 
a precision of O(10~ 5 ), the SUSY contributions to C^ik = 3,4,5,6) at the scale m w are 
the same for both b — > s and b —>■ d transitions. The contributions from the gauge boson 
W ± , the charged Higgs and the charginos are 



Cf{m w ) 
Cf{m w ) 



{0.00155, -0.00197, 0.00066, -0.00197}, 
{-0.00001, 0.00004, -0.00001, 0.00004}, 
{0,0.00003,0,0.00003}. 



For Cjy(Mw) and Cag(Mw), the NLO level numerical results are 

' -0.2175 -0.0422 -0.0007 - 0.0002/ = -0.2604 - 0.0002/, 



C 7l (m w ) 



-0.2175 -0.0422 -0.0009 - 0.0002/ 



G% {m w ) 



(40) 
(41) 
(42) 



d 



-0.2606 - 0.0002/, b^s 



(43) 



C 8g (m w ) 



C 0T178 C 0_0^73 C 0_0002 

C% M (m w ) Cjf (m w ) Cgf (m w ) 

-0.1178 -0.0473 -0.0002 - 0.0001/ 

-v ' 



-0.1653, b^d 



-0.1653 - 0.0001/, b^s 



Cl g M (m w ) C jf (m w ) 



C?„ (mw) 



(44) 



The new physics contributions to Ck{M\y) (k = 3,4,5,6) are clearly two orders smaller 
than their SM counterparts and therefore can be neglected safely For Cj y (Mw) and 
Csg(Mw) the charged Higgs contribution is dominant over the chargino contribution, but 
still much smaller than their SM counterparts. Obviously the case A is not phenomeno- 
logically interesting, since the SUSY effect is too small to be separated from the SM 
contribution though experimental measurements. 
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2. Case B 



Now we turn to Case B. For this case, the SUSY contributions to C^ik = 3, 4, 5, 6) are 
still negligibly small: (a) the charged Higgs contributions are at the O(10~ 7 ) level; and 
(b) the chargino contributions are at the O(l0~ 5 ) for both b — > s and b — > d transitions. 

For C 7l (M w ) and C 8g (M w ), however, the SUSY contributions are significant: 



C 7l (m w ) 



-0.2175 -0.1128 +1.0111 + 0.0063/ = 0.6808 + 0.0063/, b -»• d 



-0.2175-0.1128 +1.0193 + 0.0091/ 



0.6890 + 0.0091/, b 



(45) 



C 8g (m w ) 



' -0.1178-0.1103 +0.4622 + 0.0007/ 

ci g M (m w ) Cjf(m m ) 



C$ g (m w ) 



^1^^1103+0.4631 + 0.0010/ 



0.2341 + 0.0007/, b 
0.2350 + 0.0010/, b 



^8g 



At the lower scale mj, they are 

-0.3067+0.5896 + 0.0039/ 



C 7l (m b ) 



-0.3067+0.5947 + 0.0058/ 



C?* (m b )+C^ (m b ) 



C 8g (m b ) 



-0.1500 +0.2449 + 0.0005/ 

Ci g M (,m b ) Cg ± (m 6 )+<7f ff ± (m i) ) 

-0.1500 +0.2455 + 0.0007/ 



0.2829 + 0.0039/, b 

0.2880 + 0.0058/, b 

0.0949 + 0.0005/, b 

0.0955 + 0.0007/, b 



From the numerical values in Eqs. 



, one can see that 



d 



(46) 



(47) 



(48 



At the scale Mw, the charged Higgs contributions to both CV 7 and C 8g have the 
same sign with their SM counterparts, and are comparable in size with them. The 
chargino contributions, however, have an opposite sign with C^' 1 and C|^ M , and 
much larger in size than them. 

At both energy scales mw and m b , the net new physics contributions to C 7l and 
C 8g are always positive and consequently cancel their SM counterpart. The total 
value of these two coefficients therefore become positive after the combination of 
the SM and the new physics parts. 

It is easy to understand why the new physics contributions in case B are much larger 
than those in case A. In case A, the new physics contributions from both charged 
Higgs boson and charginos are negligibly small. In case B, however, we have much 
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light charged Higgs boson and charginos, which leads to large SUSY contributions. 
After the cancellation among the SM and SUSY contributions, the net value of CV 7 
and C$ g is positive. 

• In case B, though CV 7 (m&) received a large supersymmetric correction and has the 
opposite sign with its SM counterpart, it's absolute value changes a little and makes 
the theoretical prediction for the branching ratio of b — > 57 decay still consistent 
with the data. The reason is rather simple: the branching ratio Br(b — > sj) is 
basically proportional to |C7 7 (m&)| 2 . 



C. Branching ratios: data and theoretical prediction 



Using the decay amplitudes as given in Refs. |29l? ] and the coefficients a, in Eq. (JHHj) . 
it is straightforward to calculate the branching ratios of those twenty one B — > PP decay 
modes in the SM and mSUGRA model. 

In order to show more details about the ways to include the SUSY contributions, we 
present here, as an example, the calculations for the branching ratio Br(B — > tt + K~). 



The decay amplitudes of B L 



71 



K- 



arc 



A f (B° 



71 



A ann (B° -> 7T+K- 



-&f k F B - 
V2 



(m 2 B - ml) 



x {V ub V* s [01 + a u A + < + « + al)Ri] 



+V cb V c * s [a c 4 + a c 10 + K 

-i^ fMk ^-v tb v t : 



a 



8^4]} , 

1. 
-( 
2 



b 3 (K7T)^bf W (K,7T) 



(49) 
(50) 



with 



Ra 



2m 2 K 



(m h - m u ){m u + m d ) 



(51) 



annihilation contributions [2 
A ) on the "tree" contribution 



where the coefficients a\ have been given in Eq. (|HHj) . the coefficients bi(P\, P2) describe the 

Because of the strong Cabbibo suppression (|KifeKTsl 2 oc 
the a± term ), the four B — > ttK decays are QCD penguin 
dominant decay modes, and strongly depend on "large" coefficients a\ and a\. 

We follow the same mechanism as described in Refs. [23L 126^ to include the SUSY 
contributions to B — > PP decays. 

As mentioned previously, the SUSY contributions to the Wilson coefficients of the 4- 
quark penguin operators are very small and have been neglected. The large new magnetic 
penguin contributions in mSUGRA model can manifest themselves as radiative corrections 
to the Wilson coefficients 64^,8,10 ( or equivalently to a p - T with j = 4, 6, 8, 10 and p = u,c) 
and contained in the quantities Pp 2 2 , Pf>f 2 



, P ,EW p p^ p p,EW 



P 2 ,3 



For B — > tt + K decay, for example, the quantities P V K2 an d Pk^ W can be written as 
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pp 



4 m b 2 , . 
3ln- + --G fc ( Sp ) 



+ C 3 



3 /i 3 



D p,EW 
K,2 



+ (c 4 + c 6 

-6C 8 e f (1 + of + o$) 
(d + iV c C 2 



^ in !^ _ 3Gk (o) - G x ( Sc ) - G K (1) 

3 jJL 



4 m 6 2 
-In — + - - G K (s p ) 
3 /i 3 



9C 7 e f (1 + of + of) 



(52) 
(53) 



where s u = 0, s c = m 2 c /m1 are mass ratios involved in the evaluation of penguin diagrams, 
0.3 ± 0.3 and aif = 0.1 ± 0.3 are Gegenbauer moments for K meson 



of 



C> 



n 



3 Cr 



Ca and Ct„ = C$ a + C 5 are the so-called "effective 



r< e ff 

°77 



" Wilson coefficients. 

The explicit expressions of the functions Gk(0), Gk(X) and Gk(s p ) can be found easily 
in Ref . [26^ . The twist-3 quantities Pp 3 and Ppf^ receive the SUSY corrections in the 

same way as Pp 2 2 and Ppf^ ■ 

From Eqs. (j33l52l53p and the numerical results as listed in Table ITTT1 one can see that 

• After the inclusion of SUSY contributions, the effective coefficients C^i and 
changed their sign from negative to positive. The real parts of the coefficients aP 4I 
and consequently changed by about 60% and 10% respectively, but the 
imaginary parts of a? j remain unchanged. 

• Since the magnitude of coefficients a p A 7 and a p & j is larger than that of a 8 7 and a p l0 j 
by one or two orders, the new physics contributions to C$ g dominate the total new 
physics corrections. 

• Since only the coefficients a p - T for j = 4, 6, 8, 10 receive the SUSY contributions, 
one naturally expect a moderate or large new physics corrections to those penguin 
dominated B meson decays, such as B — > Ktt and B — > Ki]' processes. The tree- 
dominated decay modes, for example B — > tttt decays, remain basically unaffected. 

For the phenomenologically interesting B — > Kr( decays and other penguin dominated 
decay modes studied here, the large SUSY contributions will be included in the same way 
as for B — > n + K~ decays. 

Among twenty one B — > PP decay modes considered here, twelve of them have been 
measured so far. The individual measurements and the world average for the branching 
ratios of these decays are shown in Table IIV1 

In Table we show the theoretical predictions for the CP-averaged branching ratios 
for B — ► PP decays in both SM and the mSUGRA model (case B), assuming \i = 
m&/2, nib and 2m,b, respectively. And Br^ +a and Br? denote the branching ratios with or 
without the inclusion of annihilation contributions, respectively. It is evident that some 
decay modes have strong \x— dependence, and the annihilation contributions can also be 
significant for B — > K-n and B — ► Krj' decays. In the following subsections, we present the 
numerical results and show the dominant theoretical errors induced by the uncertainties 
of input parameters, and focus on those well measured decay channels. 
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TABLE III: The coefficients 0) , (fi) and ajj (j = 4, 6, 8, 10 and p = u, c) for B -> irK 
decays in the SM and the case B of the mSUGRA model. 





/i = m fe /2 


u = 


nib 




2m^ 




SM 


mSUGRA 


SM 


mSUGRA 


SM 


mSUGRA 




-0.276 


+0.221 


-0.270 


+0.325 


-0.258 


+0.422 


cg'G*) 


-0.155 


+0.058 


-0.142 


+0.104 


-0.130 


+0.145 


a\j x 10 3 


-24.3 - 17.4* 


-41.7-17.4* 


-23.9 - 14.4z 


-39.8 - 14.4* 


-22.6-12.3* 


-37.4-12.3* 


a%j x 10 3 


-31.4- 12.lt 


-48.8 - 12.lt 


-29.0 - 10.4* 


-45.0 - 10.4? 


-26.9 - 9.0* 


-41.7-9.0* 


a ii x 103 


-54.0 - 15.8* 


-58.3 - 15.8* 


-40.2 - 13.6* 


-43.9 - 13.6* 


-32.3-11.9* 


-35.7-11.9* 


a h x 103 


-60.4 - 2.9? 


-64.5 - 2.9* 


-44.8 - 3.8* 


-49.0 - 3.8* 


-36.0 - 4.0* 


-39.5 - 4.0* 


a h x 104 


4.7-0.6* 


3.3-0.6* 


3.0 - 1.1* 


1.4 -l.lt 


1.8 - 1.4* 


-0.1 - 1.4* 


a h x 104 


4.6-0.3* 


3.2 - 0.3* 


2.7-0.5* 


1.1-0.5* 


1.4-0.6* 


-0.5 - 0.6* 


«ioj x io 4 


-12.0 + 12.4* 


-18.0 + 12.4* 


-13.0 + 8.6* 


-19.9 + 8.6* 


-14.9 + 6.3* 


-23.0 + 6.3* 


«ioj x !0 4 


-12.3 + 12. 5* 


-18.1 + 12.5* 


-13.4 + 8.9* 


-20.3 + 8.9* 


-15.4 + 6.6* 


-23.0 + 6.6* 



D. B — ► 7T7T and Kir decays 

The three B — > tttt decays are tree-dominated decay modes. The central values and 
the major errors of the branching ratios (in units of 10~ 6 ) in the SM and mSUGRA model 
are 

1 \9.5±0.3(^)l|j(Fo)iJ;g(7) in mSUGRA, [b > 

Br(B- - vr-vr ! - I ° ^ ^o) + ^ in SM ' M\ 
1 I 6.3 t° 2 M T 9 (F ) i°;°( 7 ) in mSUGRA, ^ 



' ' " " 1 0.39 j™^) r»F ) 1^!( 7 ) in mSUGRA, 



where the three major errors are induced by the uncertainties rrib/2 < \i < 2m&, -F B 
0.28 ± 0.05 and 7 = 60° ± 20°. 

Fig. El shows the 7 dependence of the branching ratios for three B — > 7T7t decays. The 
dots and dashed curves correspond to the central values of the theoretical prediction in 
the SM and mSUGRA model 3 , respectively. The horizontal slashed bands show the data 
as given in Table ITVl 

From Fig. Eland the numerical results as given in Eqs. (j54H5T)|) . one can see that 



3 The central values of all input parameters except for the CKM angle 7 are used in this and other 
similar figures. The theoretical uncertainties are not shown in all such kinds of figures. 
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TABLE IV: Experimental data of the branching ratios for B — ► PP in unit of 10 6 , taken from 
the HFAG website 0|. For B° — > decay, the BaBar's result will be used in our analysis. 



Decay Modes 


BaBar 


Belle 


CLEO 


Average 


B~ -» 7T-K U 

— n 

B -» ir+K- 
B° -» vr i?° 


22.3 ±1.7 ±1.1 
12.8+^ ± 1.0 
17.9 ±0.9 ±0.7 

11.4 ± 1.7 ±0.8 


22.0 ± 1.9 ± 1.1 
12.0 ± l-3±oJ 
18.5 ± 1.0 ±0.7 
11.7±2.3±{J 


18-8-3.3-1.8 
-, 2 q+2.4+1.2 

lz - y -2. 2-1.1 
1 n+2.3+1.2 
l»-U_2.i_o.9 
1 9 n+4.0+1.7 
lz '°-3.3-1.4 


21.8 ±1.4 

12.511:5 

18.2 ±0.8 
11.7 ± 1.4 


is — > 7T^7T 

B- -» vr^vr 
B° vr 7r° 


4.7 ±0.6 ±0.2 
5.5±^ ± 0.6 
2.1 ±0.6 ±0.3 


4.4 ±0.6 ±0.3 
5.0 ± 1.2 ±0.5 
1.7 ±0.6 ±0.2 


A c+1.4+0.5 
4 -°-l. 2-0.4 
a fi +1.8+0.6 
4 -°-l. 6-0.7 

< 4.4 


4.6 ±0.4 
5.2 ±0.8 
1.9 ±0.5 


B~ -» i^r/ 

B- -» 
— — 

lr -» k 17 

B — > K rj 


3.4 ±0.8 ±0.2 
76.9 ± 3.5 ± 4.4 
2.9 ± 1.0 ±0.2 
60.6 ±5.6 ±4.6 


5.3±iJ ± 0.6 
76 ± 6 ± 9 
< 12 

68 ± 10±| 


9 9+2.8 

80+1° ± 7 
< 9.3 

89 ±+H ±9 


3.7 ±0.7 

77.6±tl 

2.9 ± 1.0 ±0.2 

65.2l|° 


B~ — * ir~ri 
B~ — * tt^t] 
B° - A 
B° - A' 


5.3 ± 1.0 ±0.3 

< 4.5 

< 2.5 

< 3.7 


5.4+^-0 ±0.6 

< I 


1 9+2.8 

< 12 

< 2.9 

< 5.7 


< 4.0 

< 2.5 

< 3.7 


ij — > 7777 

B rjrj 
B — > 77 77 


< 2.8 

< 4.6 

< 10 




< 18 

< 27 

< 47 


< 2.8 

< 4.6 

< 10 


B~ -» 

B° K + K~ 


< 2.5 

< 1.8 

< 0.6 


< 3.3 

< 1.5 

< 0.7 


< 3.3 

< 3.3 

< 0.8 


< 2.4 

< 1.5 

< 0.6 



• For these tree-dominated decays, the SUSY corrections considered here are very 
small. 

• The theoretical predictions strongly depend on the value of the form factor Fq^. 

• For B° — > 7r°7r° decay, the theoretical prediction in QCD factorization is about five 
times smaller than the measured value and cannot become consistent with the data 
within the whole parameter space. 

• The central value of Br(B° — > 7r + 7r _ ) is much larger than the data, but can become 
consistent with the data if one uses a smaller form factor Fq~™ or a large angle 
7 ~ 120°. But a small Fq^ is disfavored by the large measured decay rates for 
B — > 7r°7r° and Kit decay modes, while a large 7 around 120° is also in conflict with 
the global fit result 40° < 7 < 78° at 95%C.L.[46| and the latest direct experimental 
measurement 7 = 81° ± 19°{stat.) ± 13° (sys.) ± 11° (model) 0. 

In the SM, the four B — > Kir decays are dominated by the b —> sg gluonic penguin 
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TABLE V: The CP-averaged branching ratios of B — > PP decays in the SM and minimal 
SUGRA model ( in unit of 1CT 6 ) by using the central values of input parameters. Brf +a and 
Br? denote the results with or without the annihilation contributions. 



Decays 


fi = m b /2 


ji = m h 


fjL = 2777,5 


SM 


mSUGRA 


SM 


mSUGRA 


SM 


mSUGRA 




Br f Brf +a 


B r f Brf +a 


B r f Brf +a 


B r f Brf +a 


Brf Brf +a 


B~ -► 7T~i? U 

B~ — > ir°K~~ 
B° -► vr+K- 

B° - 7T K° 


13.6 lO.Z 

8.1 9.3 
10.5 12.4 
4.4 5.3 


on £ no a 

zO.o zo.6 
11.8 13.4 
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diagrams, with additional contributions from b — > u tree and electroweak penguin dia- 
grams. For these decay modes, although the SM predictions can become consistent with 
the measured values after considering the still large theoretical uncertainties, but the 
central values of the SM prediction are indeed much smaller than the measured values 
even after the inclusion of annihilation contributions. In the mSUGRA model, the new 
penguin diagrams induced by new particles can contribute effectively to B — > Kir decays. 
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FIG. 2: The 7 dependence of the branching ratios of B — * tttt decays in the SM and mSUGRA 
model. The dots and solid curves show the central values of the SM and mSUGRA predictions. 
The horizontal gray bands show the corresponding experimental measurements as given in Table 
III. 



The numerical results (in unit of 10 6 ) are 
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(59) 
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where the second and fourth error are induced by the uncertainties F ( 



B- 



0.28 ±0.05, 



*K 



0.34 ± 0.05, and m s (2Gev) = (105 ± 20)Mev. 
Fig. El shows the 7 dependence of the branching ratios for four B — > Ktt decays. The 
dots and dashed curves correspond to the central values of the theoretical prediction in 
the SM and mSUGRA model, respectively. The horizontal slashed bands show the data 
as given in Table HVI 
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FIG. 3: The 7 dependence of the branching ratios of the four B — > Kir decays in the SM 
and minimal SUGRA model. The dots and solid curves show the central values of the SM 
and mSUGRA predictions. The horizontal gray bands show the corresponding experimental 
measurements as given in Table III. 



From Fig. |3] and the numerical results as given in Eqs. (j57H6T?j) . one can see that the 
SUSY contributions can provide ~ 50% enhancement to the corresponding branching 
ratios, and such enhancements can improve the consistency between the theoretical pre- 
dictions and the data effectively. The central values of the theoretical predictions for 
Br(B — > Ktc) in the mSUGRA model become well consistent with the experimental 
measurements. 

As for the ratio R n and R c as denned in section I, the SM relation R c w R n remain 
unchanged in the mSUGRA model. The central values of these two ratios are: 

R S C M = 1.17, fimSUGRA = L15j (61) 

R s n M = 1.20, R^sugra = 116 (62) 

The reason is that the SUSY contributions to the four B — > Kit decays are similar in 
nature, and thus cancelled in the ratio of the corresponding branching ratios. 
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E. B — ► Krf) decays 

The unexpectedly large branching ratios of Br(B ± — > K ± r/) and Br(B° — > K°t]') were 
reported by CLEO, BaBar and Belle Collaborations la 14 l a |9]h and have been studied in 
the SM (fj and new physics models by many authors |7l. Il4l|. 

For the branching ratios of B — > K rf and .fT ?]' decays, as can be seen from table IIVI 
and the experimental measurements are about twice that of the the central values of 
the SM predictions in QCD factorization. In the mSUGRA model, the SUSY contribu- 
tions can provide an additional ~ 30% enhancements, which play an important rule in 
interpreting the rj'K puzzle. If we also consider the effects of those dominant errors, the 
theoretical predictions (in unit of 10 -6 ) are 

Rr(R-^K-n') = / 381 -^O) -H(Fo) tlf(m s ) t\l{l) in SM, 
1 I 49.4 i&ifcx) ±»i(F ) t^(m s ) m mSl 



38.1 in SM > 



49.4 tlH in mSUGRA, 

Rr(R o ^ K o W) / 35.4 tIM tlt(Fo) ^ 7 9 (m s ) ± 0.3( 7 ) in SM, 
( * ] ~ I 45.7 t 9 7 i(Fo) t^(m s ) ± 0.3( 7 ) in mSl 



mSUGRA, 

(63) 



35.4 tl 6 7 2 in SM, 

45.7 t\l;l in mSUGRA, 



mSUGRA, 

(64) 

where the individual errors are added in quadrature. The relation between Fq~* v ) and 
have been defined in Ref. 0|. It is evident that the theoretical predictions in the 
mSUGRA model are consistent with the data within one standard deviation. 

For B~ — > K~r] and B° — > K°1] decays, the annihilation contributions are less than 
2%, while the SUSY enhancements are about 30%. The numerical results (in unit of 10~ 6 ) 
are 



' " 1 3.6 t° °M t\i(F ) tki(m s ) t° 5 r(l) in mSUGRA, 
„_ . j 2.0 i°» t° l(F ) tli(fn s ) ±g-i(T) in SM, 

" 1 2.8 ±g5(A*) tH(F ) tlttms) ±0.2( 7 ) in mSUGRA, 



The theoretical predictions in both the SM and mSUGURA model are all consistent with 
the data within one standard deviation. But the consistency between the theoretical 
predictions and the data is clearly improved by the inclusion of the SUSY contribution. 

In Fig. 15 we show the 7 dependence of the branching ratios for four B — > Kv^ ) decays. 
The dots and dashed curves correspond to the central values of the theoretical prediction 
in the SM and mSUGRA model, respectively. The horizontal slashed bands show the 
data as given in Table ITVl 

F. B — > nrf) and B — > rjTfO 

These seven decay modes are tree-dominated decay processes, and new physics en- 
hancements due to the SUSY contributions are less than 10%. For the measured 
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FIG. 4: The 7 dependence of the branching ratios Br(B — > Krj^) in the SM and minimal 
SUGRA model. The dots and solid curves show the central values of the SM and mSUGRA 
predictions. The horizontal gray bands show the measured value as given in Table III. 



B — > 7r 77 decay, the theoretical predictions in the SM and mSUGRA model are 

1 V) " ' I 4.7 + 0.10*) ^|(Fo) i°J(m s ) li;°( 7 ) in mSUGRA ^ j 

and consistent with the experimental measurements within one standard deviation. For 
B~ — > 7T — 77' decay, the theoretical prediction for its branching ratio is similar with that 
for £?~ — > tt~7] decay, and may be observed soon by the B-factory experiments. 

For B° — > 7r°7/ ) and 5° — > 7/ '7/ ' decays, the theoretical predictions in the SM and 
mSUGRA model are of order 10 -7 and smaller than the experimental upper limits. 

For B meson decays involving an 77 or 77 as at least one of the two final states, some 
specific contributions such as the color singlet contribution have been discussed in |2§]. 
These contributions will be in favor of accounting for the experimental data. However, 
large uncertainties go with them. In our calculations, such contributions are not taken 
into account. 
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G. B -> if if decays 

For the three B —>■ KK decay modes, only experimental upper limits are available 
now. The B° — > K + K~ decay receives only the weak annihilation contribution. It's 
branching ratio is strongly suppressed in the QCD factorization approach. In Ref . |48| . 
the authors calculate this decay mode by employing the PQCD approach and also found a 
small branching ratio. From this decay mode, we can obtain useful information about the 
long-distance final state interaction and soft annihilations when the precise experimental 
measurement becomes available in the future. 

For B° -> K°K and B ± -> K±K° decays, they are penguin- dominant and the SUSY 
contributions can provide ~ 20% enhancements to their branching ratios, and still within 
the experimental upper limits. 



H. Uncertainties of theoretical predictions 

From the numerical results as given in Table IVI and in Eqs. (j54H6Tj) . one can see that 
the theoretical predictions still have large uncertainties. 

For most B — > PP decays, the dominant error comes from the uncertainties of the 
corresponding form factors, since the branching ratios are generally proportional to the 
square of the related form factors. The measured Br(B° — > 7r + 7r~) = (4.6 ± 0.4) x 
10~ 6 prefers a smaller F^^fo), but the large decay rates for B° — > 7r°7r° needs a large 
F ( f~ i " 7r (0). The measured large branching ratios for B — > K-k and K rj' decays also favor 
large .^"^(O) and F^ - * 11 (0). Further reduction of the uncertainties of the form factors 
is essential for us to find the signal of new physics from the B — > PP decays. 

The large uncertainty of the light quark masses is also a major source of the theoretical 
errors. For B~ — > K ~i] , for example, a 44% enhancement can be obtained by varying 
m s {2GeV) from 105 MeV to 80 MeV. 

Thirdly, the CKM angle 7 has a wide scope and can bring large uncertainties to 
the theoretical predictions for some B — > PP decays in both the SM and the mSUGRA 
model. Of course, one can also constrain the angle 7 from the experimental measurements 
of B — > Ktc decays |49| . 

The 7 dependence of the branching ratios for those measured B — ► PP decays are 
illustrated in Figs. I2HU In these figures, the dots and dashed line shows the SM and the 
mSUGRA predictions, respectively. The theoretical uncertainties are not explicitly shown 
here. 

From Figs. |2HH one can see that some decay modes (B — > 7r + 7r~, tc°K ~, tc + K ~, etc.) 
are sensitive to the angle 7 in both SM and the minimal SUGRA model, while other 
decays such as B — > tt~K °, tt°K and K rj have a weak dependence on the angle 7. By 
analyzing the expressions of the decay amplitudes, we find that if the term proportional 
to V u b is dominant over other terms in the total decay amplitude of a given decay, the 
branching ratio of this decay will has a strong dependence on the angle 7. 

Fourthly, the endpoint divergence of X H in the hard spectator scattering can produce 
large uncertainty to the theoretical calculations. But it is generally not important for 
those tree- or penguin-dominated decay processes because of the strong suppression of 
the a s and N c . 
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In the QCD factorization approach, the annihilation contributions can not be calcu- 
lated reliably, but estimated with large uncertainty. For B — > PP decays, the annihilation 
contribution may be strongly power suppressed as discussed by Ali et al.j23[. Of course, 
such assumption has given rise to some controversy. 

Finally, when considering the branching ratios in the minimal SUGRA model, 
different numerical results can be obtained by varying the SUSY parameters 
(m , mi/2, tan j3, A Q , Sign^fi)) around the given values as listed in Table EJ In this paper, 
we considered two typical sets of SUSY parameters which are still allowed by the data of 
B — > X<,7 and other measurements. In case A the SUSY contribution is small and can 
hardly change the SM predictions. On the contrary, in case B the SUSY contribution 
is significant in size and provides favorable enhancements to the branching ratios of the 
penguin dominant decay modes. 

V. SUMMARY 

In this paper, we calculated the SUSY contributions to the branching ratios of B — > PP 
decays in the framework of the mSUGRA model by employing the QCD factorization 
approach. 

In Sec. HH a brief review about the mSUGRA model was given. In Sec. IHIl we evaluated 
analytically the new penguin diagrams induced by new particles ( gluinos, charged-Higgs 
bosons, charginos and neutalinos), and obtained the analytical expressions of the SUSY 
contributions to the Wilson coefficients. The calculation of B — > PP decays in the QCD 
factorization approach is also discussed in this section. For the mSUGRA model with 
the mixing matrix as given in Eq. (J28j) . we found that (a) the SUSY corrections to the 
Wilson coefficients Ck {k = 3 — 6) are very small and can be neglected safely; (b) the 
leading order SUSY contributions to the Wilson coefficients C-j^Mw) and Cs a (Mw) can 
be rather large, and even change the sign of the corresponding coefficients in the SM. 

In Sec. IIV1 we calculated the branching ratios for twenty one B — > PP decays in the SM 
and the mSUGRA model, and made phenomenological analysis for some well measured 
decay modes. From the numerical results, we find following general features about the 
new physics effects on the exclusive charmless hadronic B — > PP decays studied in this 
paper: 

1. For those tree-dominated decays, such as B — > nn, the possible SUSY contributions 
in mSUGRA model are very small and can be neglected safely. 

2. For those penguin-dominated decay modes, the SUSY contributions to their branch- 
ing ratios can be significant, around 30 — 50%. 

3. For the four B — > Kir decays, the SUSY contributions to the branching ratios play 
an important rule to improve the consistency of the theoretical predictions with the 
data. 

4. For B — > Krj' decays, the theoretical predictions for branching ratios become con- 
sistent with the measured values within one standard deviation after the inclusion 
of the large SUSY contributions in the mSUGRA model. This is a possible inter- 
pretation for the so-called Krf puzzle. 
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5. The theoretical predictions in both the SM and mSUGRA model still have large 
theoretical uncertainties. The dominant errors are induced by the uncertainties of 
the form factors F^ p , strange quark mass m s , the low-energy scale /i ~ m h and 
the CKM angle 7. 
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APPENDIX A: ONE-LOOP FUNCTION AND THE COUPLING CON- 
STANTS 

In this Appendix, the explicit expressions of fi(x) functions and the coupling con- 
stants appeared in Eq. (fT^|) -(j2^ j) are presented. For more details, one can see Ref. 35|| and 
references therein. 



fUx) = — — — - (x 3 - 6x 2 + 3x + 2 + 6xlnx) , (Al) 
12(x — l) 4 ' 

f 2 {x) = 1 (2a; 3 + 3a; 2 - 62 + 1 - 6a; 2 In x) , (A2) 

h(x) = 2{ J_ 1)3 {x 2 - 4x + 3 + 21nx) , (A3) 

U{x) = 2{x l _ if (a; 2 - 1 - 2* In a;) , (A4) 

f 5 (x) = ——^ — -r(7a; 3 - 36x 2 + 45a;-16-121na; + 18a;lnx), (A5) 
36(a; — l) 4 

f 6 (x) = —— — - (37- 171a; + 207x 2 - 73a; 3 + 3 In a; - 81x 2 lna; + 54a; 3 lnx) ,(A6) 
54(x — l) 4 

f 7 (x) = —7— — - (-llx 3 + 18x 2 -9x + 2 + 6x 3 lnx) . (A7) 
18 — 1) 
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For the coupling constants Ti 



G(LR)^C(LR)'^N{LR)i We naVe 



( r G/?)/ 



- (r D )f 3 , 



V^^iy sin /3_ 



ml 



\/2m w cos /? 



K. 



kj j 



1 

i 

7! 



m 



mw cos (3 



cos /3 



(A8) 
(A9) 

(A10) 
(All) 

, (A12) 
(A13) 



where the K is the CKM matrix. 
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